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' Abstract. I review recent progress in perturbative QCD on two fronts: extending next-to-next-to- 

' leading order QCD corrections to a broader range of collider processes, and applying twistor-space 

(N ■ methods (and related spinoffs) to computations of multi-parton scattering amplitudes. 
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A method is more important than a discovery, 
since the right method will lead to new and even more important discoveries. 

- L. D. Landau 



INTRODUCTION 



Asymptotic freedom [ 1], for which Gross, Politzer and Wilczek received the 2004 Nobel 
Prize in Physics, provides the conceptual framework for applying perturbative QCD 
to short-distance-dominated problems in hadronic physics, such as the deep-inelastic 
(DIS) scattering process, the focus of this series of workshops. Supplemented with the 
notion of factorization |2f], and the experimental determination of parton distributions, 
^ ■ perturbative QCD has become the basis for all quantitative theoretical predictions for 

large-transverse momentum processes in hadron-hadron and ep collisions, as well as jet 
k> ; production in e+e annihilation. 

■ One might have thought that by now, with the aid of computers, perturbative QCD 

^ ' should have been "reduced to quadratures", that is, to a simple exercise in tabulating and 

numerically evaluating Feynman diagrams. Yet it is often the case that the experimental 
precision exceeds the theoretical uncertainties, due to unknown higher-order terms in the 
perturbation series. Here I will cover two topics in computational perturbative QCD, for 
which there has been a great deal of progress, although much still remains to be done. 

The first topic concerns next-to-next-to-leading order (NNLO) corrections to collider 
processes, also the subject of another talk at this workshop [3]. NNLO computations 
have been available for a limited number of collider observables for many years, but 
only now are the prospects becoming good for extending them to a broader range of 
important precision processes at hadron colliders. 



^ Talk presented at the 13th International Workshop on Deep Inelastic Scattering (DIS 05), Madison, 
Wisconsin, April, 2005. Research supported by the US Department of Energy under contract DE-AC02- 
76SF00515. 



The second topic is a rapidly developing one, in which insights gleaned from the 
topological string in twistor space proposed by Witten [4], and further developments, 
promise to provide efficient means for computing tree-level and one-loop QCD ampli- 
tudes with a large number of external partons, as well as vector and Higgs bosons. These 
amplitudes are needed for next-to-leading order corrections to a variety of processes. 



PROGRESS AT NNLO 

For most observables, the QCD perturbation series is a slowly converging one. (Techni- 
cally it is an asymptotic series, but rarely are there enough terms available in the series 
for the distinction to matter quantitatively.) Typical next-to-leading order (NLO) cor- 
rections for collider processes range from 20% to 100%. Clearly any kind of precision 
measurement, say at the few percent level, will require the NNLO terms in the series as 
well. Examples where this precision is desirable include the determination of 

• tts via jet production and event shapes in e+e^ annihilation (as well as in ep 
collisions) 

• parton distributions via DIS, Drell-Yan production, and high-/)/- jet production at 
hadron colliders 

• electroweak parameters, such as Mw, via W and Z production at hadron colliders 

• the "partonic luminosity" at the LHC 130 

• Higgs couplings. 

The progress of NNLO computations for collider processes can be charted in terms 
of the number of physical scales present in the parton-level cross sections. The more 
scales, the more difficult the computation, but the more flexible the applications. In per- 
turbative QCD with massless quarks, all relevant scales are associated with the external 
kinematics. (The dependence of the partonic cross sections on the renormalization and 
factorization scales can be determined with relatively little effort, so it can be neglected 
in this counting.) Also, dimensional analysis can be used to remove an overall dimen- 
sionful scale from the problem, leaving just the number of dimensionless ratios. 



No-scale problems 

For example, in the total cross section for e~^e~ annihilation into hadrons, or equiva- 
lently the ratio Rg+g-{s) = o{e^e^ hadrons)/a(e+e^ /i+/i^), the only physical 
scale is s, the square of the center-of-mass energy. This scale can be removed trivially, 
so Rg+e (s) is really a no-scale problem, from the computational point of view. That is, 
each term in the perturbative series for Rg+e is (for fixed renormalization scale ji) a 
pure number. Related to the lack of other scales in the problem is the totally inclusive 
nature of the observable; that is, it sums over all hadronic final states with no constraints. 
This sum can be performed using unitarity, or the optical theorem, as illustrated in fig. [H 
transforming the problem into the computation of the imaginary part of the virtual pho- 
ton propagator, or two-point function. 




No-scale processes were the first to be computed at NNLO in perturbative QCD, in tlie 
early 1990s. Besides Rg+e and the closely related problem of the semi-hadronic width of 
the T lepton, r(T ^ + hadrons) various DIS sum rules were also evaluated at this 
order: the Bjorken sum rule for neutrino scattering, /J dx [F^^{x, Q^) —F^^{x, Q^)] 0], 
the Bjorken sum rule for polarized electroproduction, and the Gross -Llewellyn Smith 
sum rule for neutrino scattering [8]. The integrals over x not only remove dependence 
of the observables on parton distribution functions, but they reduce the computation to a 
no-scale, propagator-type problem very similar to Re+e ■ 

The technology underlying all these computations was integration by parts (IBP) J^. 
Total derivatives of multi-loop integrals in D = 4 — 2£ space-time dimensions (i.e., using 
dimensional regularization), such as 



d^pd^q ...^ -y— 

dqt^ p^q^{p + q)^... 



J d^pd^q. 



2k -q 2k-{p + q) 



p^[q^]^{p + qY... p^q^iip + q)^]^ ■ ■ ■ 



(1) 



where p and q are loop momenta, and k is an external momenta, can be re-expressed as 
linear equations relating loop integrals with propagators raised to different powers. The 
absorptive parts of four- loop propagators can be related to three-loop propagators via the 
R* operation [10]. For the three-loop propagator problem, a recursive solution to the sys- 
tem of linear IBP equations, in terms of a small set of irreducible, "master" integrals, was 
implemented in the program MINCER lITlll . making possible the previously-mentioned 
NNLO results. 

The pure numbers encountered at NNLO in no- scale problems have a very simple 
analytic structure. The only algebraic quantities appearing, besides rational numbers, are 
the Riemann zeta values, C,{n), for n < 5. The results from the early 1990s stood as the 
computational state-of-the-art for many years. (Very recently, a method for handling the 
four- loop propagators with all possible topologies has proven successful [12], suggesting 
that the N^LO results for R^+g may be available before long.) They have led to (X.s{Q^) 
determinations with the smallest theoretical uncertainty — if can be made large 
enough experimentally, for example, at the iP pole. However, the experimental precision 

is highly stressed by the leading "1" in R^+g 1 + ttj/^rH : A 3% measurement 

of as{Mz) ~ 0.120 requires a parts per mil measurement at the Z pole of r(Z 
hadrons)/r(Z ji^fi^) [.13.1 . Observables beginning at order tts, e.g., e^e^ event- 



shape variables, are less demanding in this way, motivating their NNLO computation, 
which is a multi-scale problem, only now approaching completion. 



One-scale problems 

Also around 1990, the first NNLO computation of a one-scale collider process was 
carried out, the total cross section for inclusive production in hadronic collisions of a 
lepton pair via the Drell-Yan process, i.e. via a vector boson V = y, W or Z lil4il . At 
the parton level, the process pp ^ V +X, where V has mass My, introduces the single 
dimensionless ratio z = My/s, where s is the squared partonic center-of-mass energy. 
Whereas the NLO correction to the total cross section was sizable, at NNLO, for W or 
Z production at the Tevatron or LHC, the perturbative series nicely stabilized. 

The NNLO Drell-Yan result was followed quickly by the Wilson coefficient functions 
C,(z) for DIS structure functions [15] — except for the longitudinal structure function 
Fl, which begins at one order higher in a^, and whose computation was just completed 
this spring lil6l] . 

In the past few years, the NNLO corrections to two additional one-scale collider 
processes were attacked. First, the total cross section for inclusive production of a Higgs 
boson in hadronic collisions, pp ^ H + X was computed in the large nij approximation. 
In this limit, Higgs production is kinematically very similar to the Drell-Yan process, 
because V and H are both massive color-singlet particles, and no other mass scales 
remain in the problem — other than the overall Higgs coupling strength, dictated by 
the operator C{mt)Htr{G^vG^^). Indeed the first Higgs computation, via a high-order 
expansion in 1 — z (where now z = Mjj / s) ifnll was also applied to the Drell-Yan case, 
revealing a numerically small correction to the original results. 

A second Higgs production computation [18] exploited unitarity to express the par- 
tonic Higgs cross section as a forward scattering process, as shown in fig. El In this case, 
the state that scatters forward is not a single massive virtual photon, but a pair of mass- 
less initial partons. Also, not every cut is considered, but only those that cut through 
the Higgs particle (or vector boson V). The advantage of this approach is that the large 
number of phase- space integrals that have to be performed (only one example of which 
is shown on the left-hand side of fig. El can be traded for multi-loop integrals, to which 
the IBP method can be applied in an automated fashion [1^, in order to reduce the inte- 
grals to a manageable set of master integrals. In contrast to the no-scale examples, now 
all the master integrals depend on z. They also depend on the dimensional regularization 
parameter £, and have to be expanded in a Laurent expansion around e = 0, beginning 
at order 1 /e^, due to infrared divergences in the integrals. Fortunately, the IBP method 
also provides a way to determine the z-dependence of each coefficient in the Laurent 
expansion: Taking a derivative with respect to z produces an integral which can also be 
reduced to master integrals, thus generating a coupled set of differential equations lEoll 
which are readily solved in terms of special functions. In this way the exact dependence 
of the NNLO partonic Higgs cross section on z was determined [18] (see also ref. ^2)\\ ). 





FIGURE 2. Unitarity relates the gg +X cross section to forward-scattering loop integrals. 



For the Drell-Yan or Higgs total cross section, the special functions that appear are 
poly logarithms of the form Li„ (z) , defined by Li i (z) = — In ( 1 — z) , and for n > 1 , 



(2) 



In the Drell-Yan/Higgs application, one needs n < 3; the argument z may be replaced 
by a few other rational functions of z. The existence of the single scale z allows the 
analytical complexity, in terms of the number of different algebraic objects present, to 
grow significantly with respect to Re+e » but it is not yet out of hand. 

In principle, no NNLO computation of a collider process is complete without an evo- 
lution of the parton distributions from low to high scales at the same NNLO accuracy. 
Last year saw the long-awaited completion of the NNLO corrections to the DGLAP evo- 
lution kernels for parton distribution functions sl'ZR .Pijix) . These were also computed by 
considering a forward scattering problem, with a virtual photon and a quark in the initial 
and final state. (For the gluon evolution kernel Pgg{x), and also for Pqg{x), a fictitious 
Higgs-like scalar ^, coupling to gluons via ^tr(G^vC^^), is used instead of a virtual 
photon.) After renormalization and subtraction of coUinear divergences from lower loop 
orders, the NNLO evolution kernels are identified from the remaining l/e poles in the 
expression, which must be subtracted in the MS scheme for defining parton distribution 
functions (or equivalently, leading twist operators). As a by-product of the computation, 
the finite, order e*^ terms give the N-^LO contributions to the DIS structure function F2 
and the NNLO contribution to Fl h(i\ . 

For these results, somewhat more complicated special functions are required, such as 
Li„(.x;) with n = 4 for Pij{x) and n = 5 for F2 and F^. However, the ordinary poly log- 
arithms are not sufficient; a suitable generalization, harmonic poly logarithms ll23ll . can 
be used instead. In fact, these computations were not performed in terms of the variable 
X, but rather the variable appearing in the Mellin transform, f{N) = /q dx x^^^ fix)- 
In A^-space, the special functions that appear at NNLO are harmonic sums, such as the 
one-dimensional sum Sk{N) = Y!iLi and the multi-dimensional generalization of it, 



}m,mi,m2,...,mk 



(N) 



N 
i=l 



'mi,m2,...,mi 



(0 



I" 



(3) 



Although it might appear that in Mellin-moment space a no- scale problem has been 
recovered, it is of course an infinite set of no-scale problems. Indeed, at fixed moment 



A'^, the program MINCER can be used to compute anomalous dimensions at NNLO for 
and even up to A'^ = 16 [25]. However, for the case of arbitrary A'^, 
new integral reduction algorithms had to be developed [22]. 

Recently, the Mellin moments of the NNLO Drell-Yan and Higgs production cross 
sections were presented in terms of harmonic sums of the form Q [|26ll . putting the 
mathematical structure of the two types of NNLO one- scale problems discussed here on 
the same footing. 

There have been important phenomenological applications of these one-scale results. 
Also taking into account the resummation of large threshold logarithms from multi- 
ple soft-gluon resummation, through next-to-next-to-leading logarithmic accuracy iEtIi . 
the uncertainty on the total Higgs production cross section at the Tevatron and LHC 
has been considerably reduced, from perhaps 30-40% at NLO, to perhaps 10-20% at 
NNLO-i-NNLL. Some of the immediate phenomenological impact of the NNLO evo- 
lution kernels and coefficient functions on "DIS-driven" MRST parton fits ll28h was 
discussed elsewhere at this meeting 1.3- .29.1. 



Two or more scale problems 

Although the total cross sections for inclusive Drell-Yan and Higgs production are 
now known relatively well theoretically, in practice such quantities are not measurable 
experimentally. Numerous experimental cuts must be imposed to extract a signal from 
the background, for example, cuts on the transverse momentum, rapidity, and isolation 
of leptons or photons visible in the final state. In an ideal world, a flexible hadron- 
level Monte Carlo program, accurate to NNLO, would allow the effects of such cuts to 
be assessed. However, even fixed-order computations of generic ("multi-scale") NNLO 
observables are not yet available. 

A special case which can be handled in the style of the previous section is the 
distribution in rapidity Yy of a Drell-Yan pair or a Higgs boson. The former process 
can be used to extract parton distribution information from fixed-target production, or 
monitor the "partonic luminosity" via W and Z production at the LHC [iSJ, because 
at leading order it is proportional to q{xi)q{x2) with xi^2 = {My / y/s)e^^^ . Compared 
with the total cross section approach of ref. [18], a 5 function needs to be inserted into 
the phase-space integration, of the form 5{Y — Yy). The IBP method still works llsoll . 
reducing the phase-space integrals to a set of master integrals depending on z and u = 
{x\ / X2)e~^^^ , which can again be determined by integrating differential equations. The 
NNLO results have much lengthier expressions than the one-scale answers. They involve 
polylogarithms with arguments which can be irrational functions of z and m, for example 
U2[{u - 1 - i^{Au^-z{\+uY)/z/{2u)\ . The same stability of the perturbative series 
seen for the total W and Z production cross section at NNLO, holds also bin-by-bin in 
rapidity [30]. 

For problems with more than two scales, for example event shapes and generic 
hadron-coUider processes, a flexible, fully numerical approach seems necessary. The 
major bottleneck at present comes in integrating contributions containing the emission 
of two extra partons, which have quite complicated singularities as momenta become 



soft and/or coUinear. There has been important recent progress in this direction. For 
lack of space, and because these developments were described in another talk at this 
workshop, I refer the reader to that report [,3J. 



TWISTOR SPINOFFS 

There has been a great deal of recent interest in novel methods for evaluating QCD tree 
and loop amplitudes, stimulated by Witten's topological string in twistor space |4]. In 
general, it is possible to find compact representations for amplitudes for many external 
particles, and more efficient techniques for computing the amplitudes, by making full 
use of their analytic structure (which is sometimes hidden). Twistor space iIbUi is a kind 
of Fourier transform of the usual momentum-space representation of amplitudes. Very 
often, a Fourier transform can expose simplicity. Consider the time-dependence of the 
electric field E{t) associated with light arriving from the Sun. It has a pretty random 
appearance. However, transforming to energy variables, E{t) E{(o) = J dte^'^'E{t), 
reveals spectral lines, from which the presence of helium in the Sun could be deduced. 

The twistor transform is very well-suited for describing the scattering of massless 
particles. Traditional scattering variables are the four-momentum vectors — which 
are null vectors in the massless case, kf = — and their Lorentz-invariant products, 
Sij = {ki + kjY = 2ki ■ kj. We can trade the kf for spinor variables, the right- and left- 
handed, or + and — chirality, solutions to the Dirac equation, u±{ki). A shorthand 
notation for the two-component (Weyl) versions of these spinors is, 

{Xi) a = u+ {ki) , (li) a=u- (kj) . (4) 

The trade is possible thanks to the form of the positive-energy projector for massless 
spinors, u{k)u{k) = k, or in two-component notation, 

kf{o^)aa = {h)aa = {^i)a(^i)a ■ (5) 

Instead of Lorentz-invariant products, the natural variables for massless scattering are 
spinor inner-products El, defined by 

{jl) = e''l^{Xj)a{Xi)p=u-{kj)u+{ki), [jl]=e''^{lj)M)p=u+{kj)u4ki). 

(6) 

These products are the square roots of the Lorentz products, up to a phase ^, 

{jl) = ^le'^n , [j /] = ± v/S-ye-'^^' . (7) 

The utility of these variables was recognized already in the 1980s. For example, the 
Parke-Taylor tree amplitudes IsbL l34ll are for the scattering of two negative-helicity 
gluons, labelled j and /, and n — 2 positive-helicity gluons. They are termed "maximally 
helicity-violating" (MHV) amplitudes, because tree amplitudes with fewer (zero or one) 
negative-helicity gluons vanish. In terms of spinorial variables, they take a remarkably 



simple form for any n, 

xMHVj/ ^ 4tree^i+ 2+ /- „ + \-j^^J]l CS^ 

depending only on the positive-helicity spinors A,-, not the A,. 



Twistor space and MHV rules 

The twistor transform is a Fourier transform of the A/, leaving the A/ alone. The four 
coordinates of twistor space, for each of the n particles, are (Ai,A2, /i^/i^), where /i" 
is defined by 



In order to transform the MHV amplitudes dH), following ref. ||4|] we first must multiply 
them by the momentum-conserving 5-function, which can be written, using eq. dSJ) as, 

5 = j j4^expK«(A,)a(A,)«] . (10) 

Then the transformed amplitudes are 
Af^^'^''(A,-,Ai/) = / Yldli^Mmh) [ d^xA{?ii)cxp[ixXi~Xi\ oc Y[5{^i+xXi). (11) 

The product of all the linear d functions means that the amplitude is supported on a line 
in twistor space, as shown in fig.|3ta). 

Investigation of amplitudes with three and four negative helicities (NMHV and NN- 
MHV amplitudes) revealed the pattern of intersecting lines in fig. |3tb) and fig.|3tc). It 
also led to a set of "MHV rules" for QCD tree amplitudes, which are simpler than Feyn- 
man rules fss"]. Each line in fig. |3l corresponds to an "MHV vertex", which is a clever 
off-shell continuation of the MHV amplitude dH), labelled with two negative helicities 
and the rest positive. Many Feynman vertices can be lumped effectively into a single 
MHV vertex. The MHV vertices are joined with scalar propagators, that is, factors of 
1 / p^, so that no messy contractions of Lorentz indices have to be performed. An exam- 
ple of an MHV-rules diagram, corresponding to the twistor-space structure in fig.|3tc), 
is given in fig.lH 

The efficiency of the MHV rules for gluonic amplitudes motivated their quick ex- 
tension to amplitudes with massless external fermions ['s^, Higgs bosons coupling to 
gluons via Htr{G^vG^^) in the large nit limit iIbtIi . and vector bosons (y*, W^, Z), includ- 
ing DIS multi-jet processes [38]. 

In a parallel development, the twistor structure of loop amplitudes was explored and 
exploited [13911 . The simplest structure to explain is for gluonic loop amplitudes in a 
computational "toy model" for QCD, maximally = 4) supersymmetric Yang-Mills 



(a) MHV (b) NMHV (c) NNMHV 

FIGURE 3. Tree amplitudes are supported on networks of intersecting lines in twistor space. 




FIGURE 4. Example of an MHV-rules diagram, corresponding to fig.|3lc). 



theory. One-loop amplitudes in this theory can be written as a linear combination of 
scalar box integrals; no triangles or bubbles are required. The box integrals typically 
have many legs from the amplitude clustered into a single vertex of the box. In the 
MHV case, the coefficients of the boxes are either zero, or else equal to the MHV 
tree amplitudes aJ)^^^'^' i4flll . in which case the vertices all lie on a single line in 
twistor space, as shown in fig. Ha). However, this pattern turns out to be a degenerate 
case, which is resolved in the NMHV amplitudes. Here the simplest non- vanishing 
box coefficients are those with three clusters A,5,C and a single massless leg s. Their 
coefficients are supported on three lines intersecting in a ring, as shown in fig.E^b) ll4lll . 




(a) MHV (b) NMHV 

FIGURE 5. Twistor structure of box integral coefficients for one-loop amplitudes in ^ = 4 supersym- 
metric Yang-MiUs theory. 



On-shell recursion relations 



Quite recently, another approach to tree amplitudes has been developed, on-shell 
recursion relations ill |4j,|44|]. These relations are somewhat more efficient and easier to 



generalize than the MHV rules fis*], and they have very promising implications for loop 
amplitudes as well [46]. The derivation of the relations [44] is very general, relying just 
on Cauchy's theorem for functions of a single complex variable, and the factorization 
properties of amplitudes. The desired amplitude A„ =A,j(0) can be embedded into a 
family of amplitudes A„(z), labelled by a complex parameter z characterizing a shift in 
the momentum flowing through the amplitude. For example, the momenta of the pair 
of legs 1 and n can be shifted, while respecting overall momentum conservation and 
masslessness of the external legs, according to ki — > (z), fc„ Kiiz), where 

Hz) +Uz) = ki{z)+k„{z), kiiz) = kliz) = 0. (12) 

A momentum shift satisfying eq. (fT2t is best described using spinor variables, as 

Ai— >-Ai+zA„, Ai^Ai; A^^A^, A„^A„— zAi. (13) 

This is because a complex massless vector k^ has det(^) = k^ = 0. So the singular 2x2 
matrix k can be factored into a pair of spinors, as {k)aa = A^A^, as in eq. Q, where A' 
is no longer the conjugate spinor to A . 

As z varies over the complex plane, different intermediate states go on shell, gener- 
ating poles in z. Let Kij = ki+k2-\ \-ki. Then ^{z) = (^i,/ -|-zA„Ai)^ vanishes 

at 

z = z, = -Kl,/{n-\l,j\\-). (14) 

So long as A„ (z) — > as z — oo, the contour integral over a large circle C, f^-. dzAn (z) / z 
vanishes. The residue at z = 0, which is the desired amplitude A„(0), is the negative of 
the sum of the residues at z = zi- Those residues are given by the factorization of the 
amplitude into two lower-point amplitudes, evaluated in shifted, on-shell kinematics. 



The resulting recursion relation Il43ll includes a sum over /, and over the possible 
intermedate helicities h. 



A„(l,2,...,n)= £%+i(i,2,...,Z,-^i;f)-^A„_/+i(^f,,/ + l,...,n-l,n). 

h^±l=2 ' ^ij 

(15) 

These relations lead quickly to very compact forms for tree amplitudes. For ex- 
ample, there are 220 Feynman diagrams for the six-gluon amplitude. Using color 
algebra and symmetries, the information in these diagrams is represented by the MHV 
amplitudes dSj), plus two more helicity amplitudes, A6(l + , 2+, 3+,4^, 5^, 6^) and 
A6(l^,2+,3^,4+,5^,6^). Computing the first of them using the shift ([T3t . yields the 
set of diagrams shown in fig. |6l Diagram (b) vanishes because A4(— ,+,+,+) = 0. 
Diagram (c) is related to diagram (a) by the symmetry (1 ^ 6,2 ^ 5,3 ^ 4) (plus 
spinor conjugation). Diagram (a) can be evaluated in a few steps, using the MHV 





4 ■ — — T 

■ w ''^ T 

(a) (b) (c) 

FIGURE 6. On-shell recursive diagrams forA6(l + ,2+,3+,4-,5",6"). 



amplitudes, to give a single-term expression. Adding diagram (c) gives, 
-/A6(l + ,2+,3+,4-,5-,6-) 



(6-|(l+2)|3-)3 



(61)(12) [34] [45]56i2(2-|(6 + l)|5-) 

I (4-|(5 + 6)|l-)^ 

"^(23)(34)[56][61]556i(2-|(6 + l)|5-)-^ ^ 



The combination (2^1(6+ 1)|5^) = (26) [65] + (21) [15] in the denominator leads to 
an unphysical singularity in the first term of eq. ([T6b when fcg + fci is a linear combi- 
nation of k2 and ^5, which is cancelled by the second term. On the other hand, all of 
the physical factorization behavior is made manifest, in contrast to Feynman-diagram 



based representations Il32ll . The amplitude A6(l + , 2+ , 3^,4+, 5^, 6^) has three indepen- 
dent recursive diagrams. Thus the six-gluon calculation is reduced from 220 Feynman 
diagrams to 4 much simpler ones. 

The generality of this approach has led to its rapid application to different tree-level 
processes, with fermions and massive particles [45], and to one loop Essentially, 
amplitudes are being built up directly from their analytic properties. Considering that 
the heyday of 5-matrix analyticity ended with the rise of a gauge theory for the strong 
interactions, QCD, these recent computational advances may herald the final revenge of 
the analytic 5-matrix. 
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